Quark matter as dark matter in modeling galactic halo 
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Considering the flat rotation curves as input and treating the matter content in the galactic halo 
region as quark matter, we have found out a background spacetime metric for the region of the 
galactic halo. We obtain fairly general conditions that ensure that gravity in the halo region is 
attractive. We also investigate the stability of circular orbits, along with a different role for quark 
matter. Bag-model quark matter meeting these conditions therefore provides a suitable model for 
dark matter. 
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I. INTRODUCTION 

The flatness of the galactic rotation curves of neutral 
hydrogen clouds in the outer regions of galaxies has led to 
the hypothesis that galaxies and even clusters of galaxies 
are pervaded by dark matter [i|-y]. To explain the ob- 
served constant velocity, it is assumed that the decrease 
in the energy density is proportional to 1/r^, where r is 
the distance from the center of the galaxy. 

In special connection to these galactic halo and dark 
matter we would like to mention some of the previous 
works in the sequel as follows. In an investigation under 
the framework of brane-world models it has been shown 
by Rahaman et al. [J] that the observed rotation curves 
result solely from non-local effects of gravitation, such as 
dark radiation and dark pressure, and does not invoke 
any exotic matter field. Nandi et al. Q have considered 
several aspects of the 4d imprint of the 5d bulk Weyl 
radiation by combined measurements of rotation curve 
and lensing effect. In another case the energy density 
attached to noncommutative geometry has been assumed 
to diffuse throughout a region which is responsible for 
roducing stable circular orbits and attractive gravity 
For some other notatble works in this line see Refs. 



m- 



A number of candidates for this dark matter has been 



put forward [101-112 
that quark matter 
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In the present paper we propose 
^ -ll5| is such a candidate, as previ- 
ously suggested by several workers |16l - [l8| . The funda- 
mental particles of quark matter do not ordinarily exist 
as free particles since they are bound together by the 
strong interaction. However, quark matter is believed to 
exist at the center of neutron stars 16], in strange stars 
:19-22] , or even as small pieces of strange matter [23, 2^ . 

Regarding the origin and survival of quark matter, it 
is believed that the Universe underwent a quark gluon 
phase transition a few microseconds after the big bang 
[2i-l23. In fact, in 1984, Witten ^ proposed that such 
a transition at a critical temperature Tc = 100—200 MeV 
could have led to the formation of quark nuggets made 
of u, d, and s quarks at a density that would be larger 
than normal nuclear matter density. Laboratory experi- 
ments (such as CERN LHC) with relativistic nuclei aim 
to recreate the conditions similar to those encountered 
before and in the early hadronisation period [23] ■ As the 
expanding Universe cools, the hot quark-gluon-plasma 
(QGP) freezes slowly into individual hadrons [SOj . 

The survival of such nugget s has been the subject 
of many investigations [3ll - [34| . According to Bhat- 
tacharyya et al. f25'|, not only does a large number of sta- 
ble quark nuggets exist in the present Universe, but quark 
matter could be a viable candidate for cosmological dark 
matter, as already noted. Quantum Chromo Dynamics 
(QCD), the most accepted theory of strong interaction 
predicts that under extreme condition a hadronic sys- 
tem can undergo a phase transition from color confined 



hadronic matter to the QGP phase. In Astrophysics, it 
is speculated that the core of neutron stars may consist 
of cold QGP [H. 

Since we are taking the flat rotation curves as input, 
the circular orbit is assumed accordingly. The problem is 
thereby reduced to usages of the MIT bag model [3a] to 
determine conditions ensuring that gravity in the halo re- 
gion is attractive. Under those conditions, quark matter 
may be considered as one of the suitable candidates for 
dark matter. This is the motivation behind the present 
work. 

The investigations are organized as follows: in Sec. 2 
we provide basic equations whereas in Sec. Ill we have 
considered the galactic rotation curves as an input. The 
solutions of the equations have been found out for the 
cases (i) non-interacting two-fluid model, and (ii) inter- 
acting two-fluid model in Section 4. In Sees. 5, 6 and 7 
we have addressed the question of an attractive gravity 
for structure formation in the galactic region, stability 
of circular orbits and dual role of quark matter, respec- 
tively. Sec. 8 is devoted for some concluding remarks on 
the model. 



II. BASIC EQUATIONS 

In this paper the metric for a spherically symmetric 
spacetime is taken to be 



ds' 



Mr)dt^+eHr)^^2^^2^^2 



(1) 



where dfi^ = dO^ + sirP'OdcfP. We are using here ge- 
ometrized units in which G = c — 1. 

The energy-momentum tensor of the two-fluid model 
is given by 

"^o = Peff = P + Pq, (2) 

T^^T^ = -p,ff = -ip + p,), (3) 

In the above two equations ^ and ^ p and p corre- 
spond to the respective energy density and pressure of 
the baryonic matter, whereas pq and Pq to the respective 
dark energy density and pressure due to quark matter. 
The left-hand sides of equations © and ([3]) are the ef- 
fective energy density and pressure, respectively, of the 
composition. 

The Einstein field equations are listed next: 
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Si^ip + Pq) = e -" 

?>Tl{p+Pq) = e-^ 
%TT{p+Pq) = — - 



since Tl = r|. 

In the MIT bag model, the quark matter equation of 
state (EoS-1) has the simple linear form 

1 



Pq 



-{pq-m, 



iJ) 



where B, the bag constant, is in units of MeV/{fm)^ 

For normal matter, we use the following equation of 
state (EoS-2): 



p = mp, 



(8) 



where < to < 1. 

Since we are assuming the pressure to be isotropic, 
therefore, the conservation equation is 



'-^ + l-'ip.ff+Peff)=0. 



III. GALACTIC ROTATION CURVES 



(9) 



The observed flat rotation curves are often considered 
as evidence for the existence of dark matter. In such 
galaxies the neutral hydrogen clouds, observed at large 
distances from the center, are treated as test particles 
moving in circular orbits due to the gravitational effects 
of the halo. To derive the tangential velocity of such 
circular orbits, we start with the line element ([T]) and 
then observe that the Lagrangian for a test particle takes 
the following form: 
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(10) 



where as usual fl^ — 9^ + sin^9(jp. The overdot denotes 
differentiation with respect to affine parameter s. 

Since the metric tensor coefficients do not depend ex- 
plicitly on i, 9, (/), or J7, the Euler-Lagrange equation 
yields directly the conserved quantities E and L, the en- 
ergy and total momentum, respectively: E = —e'^^^'i, 
Le = r^9, and L^ — r^sin^9 <j). So the square of the total 
angular momentum is L^ ~ Lg + (L^/sin^)^. It is ac- 
tually more convenient to use Eq. ([TU| , as considered by 
Bohmer, Harko and Lobo [38[. Then the total angular 
momentum becomes L = r^fl. 

With the use of the conserved quantities E and L and 
the norm of the four- velocity u^Ufj, = — 1, the geodesic 
equation becomes 



1 



-e-'^-H^ + e^('-)f2 



+ r^{9'^ + sin'^9 4,^). (11) 



(4) 


Hence we get 


(5) 


gi/(r) + A(r)^2 _,_ gi'(r) f^^ 

v 


,(6) 


The equation of motion 




f2 + V{r) = 


)n of 


now yields the potential as 



£2 



E' 



(12) 



(13) 



V{r) - -e-^W fe-'^^-^E^ " 7^ " ^) ' ^^^^ 



Following the work of Nucamendi, Salgado and Sudarsky 
|39| , as it is more convenient to use the effective potential 
Veff{r), we write Eq. ([T2l) in the form 

Recall that for the case of circular stable orbits, the 
effective potential must satisfy the following conditions: 
(1) f2 ^ 0, (2) ^ = 0, and (3) ^^ > 0. By using 



the first condition gives directly 



E^^e'^il + ^ 



from Eq. (fT5|). while the second condition yields 

r2 1 
r^ 2 
These equations can also be rewritten as 



E' 



1-\tv' 



and 



j2_ 2 ^ ^ 



l-\rv'' 



(16) 



(17) 



(18) 



(19) 



(20) 



Before considering the third condition, Veff{r)rr > 0, 
we need to obtain an expression for the tangential veloc- 
ity u"^. In terms of proper time, this is given by [40| 






e-'r^nK 



Substituting for (l"^, we get from Eq. ([T 



e-L^ 1 , 
^ — r v . 

^2^2 2 



This expression can be integrated to yield 



(21) 



(22) 



where Bq is an integration constant and I — 2(u''^)^. 
Equivalently, the line element ([T]) can be written 



ds^ 



l^j de + e^'^^Ur^+r^d^''. (23) 



As a result, the model has a well-defined Newtonian limit 
|38| . (Observe also that Ve//(r) ~ 1/f^, characteristic of 
dark matter, as a consequence of the fiat rotation curves.) 



Now from Eq. (IT6| 
6i2 
r 



Veff{r)rr = ^ " E'e-''{iy'f + E^e-'^v". (24) 



By making the substitutions for L'^, E^, and i/, the second 
derivative reduces to 
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Veff(r)rr = "J > 



(25) 



confirming the existence of stable orbits (a detailed anal- 
ysis of which will be provided in the Sec. VI). 



IV. SOLUTIONS 



In this section we determine e^^''^ in the line element 
([T]) by considering two cases, a non-interacting two-fluid 
model and an interacting two-fluid model. 



A. Non-interacting two-fluid model 

In our first case, the two fluids, normal matter and 
quark matter, do not interact. The resulting conservation 
equations are therefore independent of each other. Using 
Eqs. (O and ©, we have 

dp , / 1 + m \ 



and 



dpq 
dr 



2v\pq - B) = 0. 



(26) 

(27) 

(28) 

(29) 

where pQ and pog are integration constants. Equations 
dH), (P, and © now yield 



The solutions are 



and 



P = Por 



-l{l+m)/2m 



Pq = B + poqr ^', 



(30) 



Stt [p(l 4- 3m) + {2pg - 4B)] = 

. r AV {ly')^ „ 2v' 

This equation is linear in e^^ and leads to 



(e~^r^+') = 87r( y )r3+'[/,(H-3TO) + 2p,-4B] (31) 



by making use of Eq. ^F2\ . Incorporating Eqs. 
(Uni), we now obtain 



and 



Dr' 



+ '-8. 



Po 



(1 + 3m)r2-'(l+™)/2m 



A + l - l{l + m)/2m 
„2-2/ 2Br2- 



Zpog?- 



4-; 

where D is an integration constant. 



4 + Z 



(32) 



B. Interacting two-fluid model 

In our second case, the two fluids are assumed to interact. 
The resulting conservation equations will then take the 
following form: 



and 



dp 
dr 



dr 



1 



2m 



2v'{p,-B) 



P = Q 



-3g. 



(33) 



(34) 



The quantity Q expresses the interaction between the 
dark-energy components. We assume that there is an 
energy transfer from quark matter to normal matter. So 
a positive Q is a natural choice that ensures that the 
second law of thermodynamics is fulfilled [41]. Moreover, 
Eqs. ([55)) and (IMl) imply that the interaction term is 
proportional to r and vanishes as r — > oo. To meet these 
criteria, we assume that 



n 1 ^n Qo 



(35) 



where < n < a for some a < 1 and Qo is a constant 
of proportionality. (This assumption is made strictly for 
computational convenience, allowing us to write explicit 
solutions; so no particular physical significance should be 
attached to the range on n). The respective solutions of 
Eqs. (|33|) and ([M)) are now given by 



P = Po{in)r 



-l(l+m)/2m 



Qor 



l-n 



1 -~n + l{l + m)/2m 



(36) 



and 



Pq — B + pog(i 



„-2; 



3Qo»' 



l-ra 



l-n + 2V 



(37) 



where po(m) and pog(ii 
The solution of Eq. 



Fr 



'2-1 



-Stt 



are integration constants. 
30|) now becomes 

Po(m)(,l + 3TO)r ^ ^ >' 



A + l-l{l + m)/2m 



2po<j(i 



^2-2; 



-«. 



2Br^ 
4:-~l 4 + 1 

Qo(l + 3m)r3-" 



[l-n + l{l + m)/2m]{5 - n + I) 



{l-n + 2l){5-n + l) 



(38) 



where F is an integration constant. Observe that apart 
from the integration constants, the first part of solution 
is the same as solution ([5^ . 



ON THE QUESTION OF AN ATTRACTIVE 
GRAVITY 



We may consider the question of an attractive gravity 
in the halo region by studying the geodesic equation 



dr2 
This equation implies that 



(h^dx^^ 
^^ dr dr 



dr2 



1 



dr 



dr 



dt\' 
d^J 



Since r = 0, as before, we get from Eq. 
d^r 1 xr. , 1 ^ f dt 



<0 



(39) 



(40) 



(41) 



as long as e~^ > 0, thereby yielding an attractive gravity 
in the halo region. 

In the non-interacting case related to Eq. (15^ . the 
dominant terms are the last two because of the size of 
the bag constant B. We can see from the EoS-1 in Eq. 
that 



1 



Pq 



(po,r-2' - 3B). 



(42) 



So whenever the pressure pq is positive, pog > 'iBr^^ > 
3B. It follows that 



^eff - 



2poqr 



2-21 



4-1 



2Br^ 
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(43) 



.-Kr) 



is positive and since po <■ Pq and / « 0.000001 ^^ 
is also positive. 

The same argument can be applied to the interacting 
case because of the similarity between Eqs. ([5^ and ([55]) . 
So we obtain an attractive gravity in both cases. 

The conclusion is somewhat different if Pq is negative. 
Based on Eq. (|42|) . poq can still be large enough to yield 
e~rf > 0. Now suppose that the corresponding expres- 
sion in Ec^. (1551) is also positive. At a first glance it seems 
difficult to quantify precisely the structure of quark mat- 
ter as much of it even now remains conjectural. However, 
for certain combinations of values of the various parame- 
ters, the extra terms involving Qo will produce a negative 
value for the total. So we no longer have an attractive 
gravity for the interacting case. Of course, if the pressure 
Pq is negative and sufficiently large in absolute value, then 
£--'*('•) < in both cases, so that quark matter would not 
be a suitable model. 



VI. ANALYSIS OF THE STABILITY OF 
CIRCULAR ORBITS 

In the Introduction we have assumed the circular orbits 
as an input of flat rotation curve and have shown stable 



via the Eq. (|25|) . Here we analyze the stabihty of the 
orbits in a more detailed form to distinguish between the 
interacting and non-interacting cases. Let a test particle 
with four velocity U" — ^^ moving in the region of 
spacetime given in (23). Assuming — 7r/2, the equation 
g^aU^'U" = -m§ yields 



(44) 







(J)"..^.m 


with 
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E^ 


(i '"''''\ i^-'(ii ^' 


[ Bo ) ' V ' r^ 



(45) 



Here the two conserved quantities, namely relativistic 
energy (E) and angular momentum (L) per unit rest mass 
of the test particle respectively are 



E = — and L = — . 
mo mo 



(46) 



If the circular orbits are defined by r = i?, then ^ — 

dV 
dr 



and, additionally, ^ \r=R= 0. Above two conditions 



result 



L = ± 



I 



2~l 



:R 



and, using L in V{R) = —E"^, we get 



2B, 



^fi''/2 



The orbits will be stable if ^ 



(47) 



(48) 



r^j^< and unstable if 



dr 



2 |r— i?-' 



""^ ' =R> 0. 

By putting the expressions for L and E in ^^ 



\r=R 



and then by using the Eq. ((38|) of the interacting case, 
finally we get 
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where 

Fi=87r(f) 

F2 = 87r(f) 
F3 = 87r(f) 

-^4 = 87r ( jj 1^ [l-n+Z(l+m)/2m](5-n+0 "" {l-n+2l){b-n+l) 

We note that for non- interacting case when Qo = 0, 
4;:^ \r=R< since, / « lO"*" and m,n < 1. Thus the 
circular orbits are always stable for non-interacting 
situation. When, Qo 7^ i.e. for interacting two fluids, 
F4 < 0. However, if Qo is very small then, ^^ \r=R 
is also negative and circular orbits are stable. The last 
condition implies that the interaction is very weak. This 
weak interaction is a characteristics of quark matter at 
higher temperatures. Thus, in the Weakly Interacting 
Massive Particle (WIMP) category of the dark-matter 
halo, besides neutrino-like light particles, quark matter 
may be one of the possible candidates, which, in turn, 
has important implications for the galactic structure 
formation and evolution. 



VII. BEHAVIOR OF QUARK MATTER: A 
DUAL ROLE 



The above results may have an interesting interpreta- 
tion on a cosmological scale. If pq and pq refer to the 
matter content in the galactic halo region, then, as we 
have seen, if Pg > 0, then bag-model quark matter be- 
haves like dark matter. Now suppose that pq refers to 
the energy density obtained by including a region be- 
yond the halo large enough so that the resulting reduced 
value leads to Pq < in Eq. (|42|) . Assume also that 
I = 2{v'^)'^ decreases so gradually in the outward radial 
direction that we may treat it as a constant. Substitut- 
ing the expression for pq from Eq. (|42p in the Friedmann 
equation [d"^ a{t) / dt^) / a{t) ~ — ^(p + 3p), we obtain 



1 d^a{t) 
a{t) dt^ 



47r, 
-^(Pg+POg'' 
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3B). (50) 



Sufficiently small pq and poq now result in a positive ac- 
celeration, in the manner of dark energy. It seems there- 
fore on the galactic level, quark matter behaves like dark 
matter and on the global level like dark energy. 



VIII. CONCLUSION 

By taking the flat rotation curves as input and treating 
the matter content in the galactic halo region as quark 
matter, we obtained a spacetime metric for the galactic 
halo. 

One can see that the solutions given in the Eq. ((32)) 
and Eq. (I55|) are the interior solution of the spacetime 
metric given by the Eq. (1231) which is neither asymptot- 
ically flat nor a spacetime due to a centrally symmetric 



black hole. The metric describes that region for which 
the tangential velocity of the test particle is constant. 

Since flat rotation curves were assumed, the problem 
of determining the suitability of quark matter as a model 
for dark matter was reduced to finding conditions under 
which gravity in the halo region is attractive. 

We used the MIT bag model whose equation of state 
is Pq = \{Pq ~ 4-B) and considered two cases, non- 
interacting and interacting two-fluid models (referring to 
quark matter and baryonic matter). In both cases, if 
the pressure Pq > 0, gravity in the halo is attractive. 
li Pq < 0, with a sufficiently large absolute value, then 
gravity in the halo is repulsive. There exist values of Pq 
between these extremes in which gravity is attractive in 
the non-interacting case but not in the interacting case. 

We have also investigated the stability of circular or- 
bits and shown that the circular orbits are always stable 
for non-interacting situation. However, interacting two 
fluids circular orbits are found to be stable under certain 
conditions which imply that the interaction is very weak. 

In our investigation it is also revealed that on the galac- 
tic level, quark matter behaves like dark matter whereas 
on the global level like dark energy. 

In summary, however, quark matter fitting the MIT 
bag model is indeed a suitable model for dark matter 



under fairly general conditions. 

As an additional remark, in the interacting case be- 
tween quark and neutral hydrogen a rise of temperature 
is expected to happen. If this environment of increased 
temperature in the interstellar medium is detectable ob- 
servationally then the present model is physically vi- 
able one. The aspect of thermodynamical processes in- 
volved in this phenomenon can be undertaken in a future 
project. 
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